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A  comparison  between  the  entirely  statistical  and  the  combined  kinetic  and  statistical 
models  for  amine-epoxy  networks,  both  studied  previously  in  the  literature  is  made.  The 
chemistry  studied  involves  the  generation  of  a  secondary  site,  which  reacts  in  the  network 
through  chainwise  addition,  rather,  than  the  stepwise  addition  mode  of  the  primary 
reaction.  It  is  shown  that  the  combined  model  ( Bokare  and  Gandhi,  1980)  is  exact  and  the 
random  model  (  Bidstrup,  1986;  Riccardi  and  Williams,  1986;  and  Tsou  and  Peppas, 
1988)  is  approximate.  The  assumptions  underlying  the  two  models  are  delineated.  It  is 
shown  that  under  certain  circumstances  the  predictions  of  the  two  models  can  differ  by 
several  hundred  percent.  Further  the  combined  model  is  extended  to  predict  the  postgel 
properties  of  the  network. 
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INTRODUCTION 


Application  of  the  classical  Flory-Stockmayer1  theory  for  polymer  growth,  and  that 
of  its  several  other  manifestations^-^  js  contingent  upon  the  reactants  meeting  certain 
conditions.  These  conditions  have  been  well  discussed^  and  ensure  that  the  reaction 
probability  between  any  two  sites  is  random.  Deviations  from  the  the  classical  paradigm, 
include,  chainwise  reactions^,  unequal  reactivity^,  substitution  effects^,  and  cyclization  in 
finite  species**.  Given  a  well  established  chemistry  it  is  possible  to  account  for  these 
non-idealities,  at  least,  to  some  degree  of  approximation^.  The  present  paper  is  concerned 
with  the  non-ideality  introduced  by  chainwise  side  reactions  in  stepwise  systems. 

The  chainwise  reaction  is  conditional  upon  the  initiation  step.  Further  non-random 
complications  appear  due  to  termination,  chain  transfer  reactions  and  diffusion 
limitations111-1 *.  "Living"  aniomc  polymerization  is  the  simplest  case  of  chainwise 
chemistry 1 1 .  Here  the  reaction  between  two  sites  is  conditional  only  upon  one  of  the  sites 
having  been  initiated.  Dusek  and  Somvarsky 1 2  have  analyzed  this  network  model 
chemistry  and  shown  the  nature  of  the  deviations. 

Consider  the  linear  polymer  case  in  homogenous  growth  conditions.  In  stepwise 
systems,  if  we  neglect  the  last  three  non-idealities  mentioned  above,  the  chain  lengths  have 
a  geometric  distribution.  In  contrast,  if  the  reaction  mechanism  is  "living"  anionic,  and  the 
initiation  step  is  instantaneous,  the  chain  lengths  have  a  Poisson  distribution1 3. 

This  work  explains  the  modelling  strategy,  where  both  the  stepwise  and  chainwise 
mechanisms  are  operative.  Specifically,  we  apply  these  results  to  the  amine-epoxy 
chemistry.  This  chemistry  has  been  well  studied  in  the  literature1 4- 16  and  several 
workers1^1  have  developed  network  models.  The  reaction  between  the  amine  hydrogen 
and  epoxy  (1)  is  usually  the  main  reaction.  However,  the  hydroxyl  unit  generated  on  the 
epoxy  can  further  react  with  another  epoxy  (2)  to  continue  the  chain  growth.  This 
secondary  reaction  is  like  "living"  anionic  polymerization1 22  por  this  reaction  the  amine 
acts  as  an  initiator  and  the  epoxy-hydroxyl  reaction  is  the  conventional  propogation  step. 
However,  in  usual  amine-epoxy  networks,  the  amines  are  stoichiometrically  in  balance 
with  the  epoxies  and  the  propogation  step  is  much  slower  than  the  initiation  step.  This 
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could  have  led  some  workers  18-20  to  ignore  the  chainwise  nature  of  the  side  reaction. 
Their  models  are  therefore  approximations  to  the  chemistry  they  assume.  Previous  work  by 
Bokare  and  Gandhi^l  accounts  for  the  chainwise  nature  of  the  side  reaction.  Dusek^  has 
also  presented  a  model  faithful  to  this  chemistry  in  the  limit  of  slow  propogation  step. 


We  shall  compare  the  modelling  assumptions  inherent  in  past  modelling  strategies. 
In  route,  we  shall  develop  a  ger  v*_  framework  for  growth  models  of  networks  with  both 
stepwise  and  chainwise  reactions.  We  shall  also  explain  the  chemical  situations  under 
which  the  deviations  due  to  chainwise  reaction  are  pronounced.  Finally  a  model  for 
calculating  post-gel  properties  for  such  networks  shall  be  presented. 

For  reasons  of  comparison  and  mathematical  ease  we  have  restricted  the  analysis  to 
amine-epoxy  networks.  Similar  considerations  hold  in  other  important  chemistries  like, 
epoxides  and  alcohols23,24>  epoxides  and  carboxylic  acids,  isocyanates  and  alcohols, 
isocyanates  and  amines. 


CHEMISTRY 


Any  model  for  network  growth  is  only  as  good  as  the  chemistry  it  assumes.  We 
shall  limit  our  model  to  equations  3-5  as  the  three  major  reactions.  Other  reactions  could 
possibly  be  important^, 26  but  they  are  not  the  subject  of  study  here.  Further,  we  shall 
ignore  the  substitution  effect  reported  on  the  diamine  (equation  4),  as  the  major  objective  of 
this  study  is  to  report  on  the  importance  of  the  third  or  the  etherification  reaction  (equation 
5),  on  the  network  models. 


primary  amine-epoxy  reaction:  STEPWISE 


(3) 


secondary  amine-epoxy  reaction:  STEPWISE  (substitution  effect  on  the  amine  is  observed; 
it  will  be  ignored  in  this  study) 


R 


epoxy  hydroxyl-epoxy  or  etherification  reaction:  "LIVING"  ANIONIC 

For  a  comprehensive  study  of  the  epoxy  chemistry  the  reader  is  referred  to  a  review 
paper  by  Rozenberg^.  On  the  importance  of  the  third  reaction,  many  workers  have 
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reported  on  epoxy-phenol  systems.  Here  the  phenol  hydroxyl  acts  as  the  initiator  for  the 
propogation  step  between  the  epoxy  hydroxyl  and  the  epoxy.  Alvey27  found  (model 
compounds  diglycidyl  ether  of  Bisphenol  A  and  Bisphenol  A)  that  even  though  the 
etherification  rate  became  slower  with  progressive  disappearance  of  epoxies,  it  was 
responsible  for  gelation  in  an  otherwise  linear  system.  Zukas,  Dunn  and  Gilbert^ 
reported  etherification  as  an  important  reaction  even  in  the  case  of  excess  amines  (model 
compounds  methoxyglycidyl  ether  of  Bisphenol  A  and  4,4-aminodiphenyl  sulfone). 
Etherification  rate  can  be  promoted  by  selective  catalysts  and  temperature^, 30  The 
phenol-epoxy  system  has  also  been  studied  by  Bantle  and  Burchard23,24  They  report 
good  agreement  with  experimental  results,  even  though  they  ignore  the  chainwise 
mechanism  in  their  model. 

For  most  amine-epoxy  systems  the  etherification  reaction  is  considerably  slower 
than  the  amine-epoxy  reaction.  Kumar  and  Gupta^2  have  suggested  a  ratio  of  around  100 
for  kj/k3  to  explain  experimental  data.  This  might  render  some  of  the  previous 
approximate  models1 8-20  in  g00d  agreement  with  experimental  studies.  However, 
situations  can  arise  in  amine-epoxy  and  similair  chemistries  where  such  approximations  are 
poor.  One  such  case  is  that  of  excess  of  epoxies  in  the  initial  stoichiometric  recipe.  Our 
working  premise  shall  assume  reactions  3-5  as  the  most  important  reactions,  without 
limiting  itself  to  any  relative  importance  among  them. 

THE  MODELS: 

The  concept  of  superspecies^  can  be  successfully  used  to  demonstrate  the 
modelling  procedure  for  this  chemistry.  Superspecies  are  defined  as  the  minimum  number 
of  distinct  units  in  the  system,  whose  combination  to  yield  the  network  can  be  statistically 
considered  random.  The  superspecies  must  be  build  kinetically.  Non-idealities  arising 
from  the  failure  of  the  assumptions  of  the  classical  Flory-Stockamyer  theory  to  hold,  can  be 
tackled  by  defining  an  appropriate  set  of  superspecies.  For  "living"  anionic  systems  the 
entire  chain  length  distribution  constitutes  the  set  of  superspecies.  Once  the  distribution  of 
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the  chains  grown  starting  from  an  initiator,  is  solved  from  the  kinetic  equations,  the 
statistical  models  can  be  used  to  combine  these  chains  together  to  form  the  network. 
Consider  a  comparison  between  the  two  models  we  wish  to  discuss  from  this  perspective. 


A.  COMBINED  MODEL:  Our  treatment  shall  use  the  model  first  presented  by  Bokare 
and  Gandhi^l.  This  model  is  exact  and  overcomes  the  approximations  made  in  some  later 
studies  17-20.  The  complete  set  of  superspecies  can  be  identified  by  considering  the 
diamine  as  an  initiator.  There  are  two  kinds  of  chains  that  are  grown  anionically, 
depending  on  the  state  of  reactivity  of  the  amine,  those  emanating  from  a  secondary  amine 
(Fig  la)  and  those  from  a  tertiary  amine  (Fig  lb).  These  combined  with  the  unreacted 
epoxy  and  unreacted  amine  constitute  the  complete  set  of  superspecies. 

1)  unreacted  epoxy  (e) 

2)  unreacted  amine  (p) 

3)  secondary  amine  sn,  with  n  epoxies  attached  on  the  side  chain  (Fig  la) 

4)  tertiary  amine  tn,  with  n  epoxies  attached  on  the  side  chain  (Fig  lb) 

( paste  Fig  1 ) 

The  upper  limit  on  n,  nmax,  measured  experimentally  is  an  indication  of  the  relative 
importance  of  etherification  rate.  Etherification  rate  k3  is  zero  if  nmax=l  for  secondary 
amine  and  nmax=2  for  tertiary  amine.  In  the  combined  model  no  upper  limit  is  imposed  on 
nmax.  Having  identified  the  superspecies,  we  follow  the  procedure  outlined  above.  We 
first  obtain  the  concentration  profile  of  the  superspecies  and  then  combine  them  statistically. 

For  complete  detail  of  the  model  and  its  mathematics  we  refer  the  reader  to  the 
original  paper^l.  For  sake  of  completness  and  since  we  shall  be  using  the  same  combined 
Bokare  and  Gandhi  model  for  predicting  the  post-gel  properties  we  briefly  reproduce  the 
mathematics  here.  The  reaction  equations  in  terms  of  the  above  based  abbreviations  for  the 
superspecies  are: 
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p  +  e 


(6) 


k3 


s" +  e~ — *■  s"*> 

n=l,2,3.... 

(7) 

Sn  +  e  — - ►  ‘n+1 

n=l,2,3... 

(8) 

k3 

n=2,3,4... 

(9) 

The  kinetic  equations  for  the  growth  of  these  superspecies  are: 

'  3e=2k‘ep 

(10) 

V(VVi>  +  k2es» 

n  >  1 

(ID 

'^e  =2k3'(v‘"-l)  '  ^es“-‘ 

n  >  2 

(12) 

where  H  is  defined  to  be  zero  and  9  is  the  time  variable. 

To  solve  this  set  of  infinite  differential  equations  we  define  probability  generating 
functions: 

oo 

s  =  2wZnsn  and  1=^2"^  (I3a&b) 

n=l  n=l 

Reducing  equations  2  and  3  in  terms  of  s  and  t  as  a  function  of  p  we  obtain: 


£ [<i+koM-*  (l4) 

dp  K0  P 

(15) 

where  K^=  and  Ko=k2/k3 

In  our  calculations  we  set  lq=k2  therefore  Ka=1/2  and  vary  Kq  These  equations 
can  be  solved  analytically  to  obtain  s  and  t  in  closed  forms.  (Equation  20  &  21)  Higher 
derivatives  can  be  used  to  obtain  Insn,  Intn,  Zn^sn  and  Xn2tn. 
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(16a-d) 


These  higher  derivatives  become  increasingly  tedious  to  compute  and  the  reader  who  needs 
them  is  referred  to  the  original  paper  21.  The  formulae  we  need  are  reproduced  in 
Appendix  A. 

We  now  use  the  recursive  technique^  to  predict  the  network  structure  and 
properties.  Looking  in  and  looking  out  of  a  randomly  chosen  epoxy,  we  write  the 
expectation  value  of  the  weight  attached  to  it  (see  Fig  1).  The  functionality  of  the  epoxy 
monomer  is  g,  therefore,  the  weight  looking  in  a  randomly  chosen  epoxy  is  the  molecular 
weight  of  the  epoxy  and  the  weight  looking  out  of  the  other  (g-1)  epoxies. 

E(We'n)  =  Mg  +  (g-l)E(Weout)  (17) 


When  looking  out  of  an  epoxy,  it  can  either  be  an  epoxy  attached  to  a  secondary  or  to  a 
tertiary  amine  or  it  could  be  unreacted: 

£(<*)  =  £  -^{(n-l)E(W;n)  +  E(wi")}  +  £^{(n-l)E(Wem)  +  E(Wam)} 
n=l  n=l 

(18) 

The  subscripts  "a"  and  ”e"  stand  for  amine  and  epoxy  respectively,  e0  are  the  number  of 
initial  epoxy  monomers.  The  first  term  covers  the  possibility  of  the  epoxy  being  any  of  the 
n  epoxies  attached  to  a  Sn  secondary  amine  and  the  second  term  where  it  is  any  of  the  n 
epoxies  attached  to  a  tn  tertiary  amine.  If  the  epoxy  is  unreacted  then  there  is  no  weight 
attached  looking  out  of  it. 

Similarly  choosing  an  amine  at  random: 
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E(Wai")  =  Ma  +  (f/2-l)E(Waout) 


(19) 


where  the  amine  monomer  has  f  hydrogens  or  f/2  amine  nitrogens,  and  p0  are  the  number 
of  initial  amine  monomers.  Also  the  expectation  weight  looking  out  of  an  amine  is  the  sum 
of  the  three  possibilities,  it  is  a  secondary  amine,  a  tertiary  amine  or  it  is  unreacted. 


E OC) 


n=l  n=l 


(20) 


Solving  for  the  expected  weight  attached  to  any  epoxy  and  amine  chosen  at  random 


E(Wa)=  (f/2)E(Wa°ut)  +  Ma  (2i) 

E(We)  =  gE(We°ut)  +  Mg  (22) 


Then  the  average  molecular  weight  of  the  reacting  system  at  any  given  conversion  is: 


fMg 


fMe  +  2pgMa 


E(Wj) 


2pgMa 


fM^  +  2pgM; 


■E(Wa) 


(23) 


where  p  =  Po/eQ  is  the  initial  recipe  ratio,  that  is,  the  number  of  amino  groups  to  the  epoxy 
groups  in  the  initial  recipe.  Ma  and  Mg  are  the  molecular  weights  of  the  monomers. 

The  condition  for  the  gel  point  is: 


l=Tl*p(g-l) 

(02-1)  < 

X  (nsn  +  ng 

n  »  1 

2 

X  (n2sn  +  n2g 

ft  *  i 

►  _< 

X  (ns» + nt»} 

n»l 

► 

l  J 

l  P^  j 

l  Peo  J 

. 

(24) 
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B.  RANDOM  MODEL:  Wc  call  any  model  that  fails  to  recognize  the  chainwise  nature 
of  the  etherification  reaction,  random.  In  these  models  18-20  the  distinct  units  that  compose 
the  network  are  identified  (Fig  2)  and  kinetic  equations  solved  for  their  concentration 
profile.  Then  the  statistical  arguments  are  employed  to  develop  the  network.  The  units 
chosen  thus,  are  not  sensitive  to  the  operative  growth  mechanism,  and  therefore  they  are 
not  random  in  their  statistics.  It  would  be  incorrect  to  call  them  superspecies,  a  name 
reserved  for  units  that  are  chemically  random  in  their  statistics.  For  bifunctional  diamines 
and  epoxies  eight  structures  can  be  isolated  as  these  units  (Fig  2).  We  label  such  models 
random,  because  they  would  have  the  correct  units  as  the  superspecies  if  the  mode  of 
reaction  was  stepwise  only. 

(paste  Fig  2 ) 

The  details  of  the  modelling  strategy  for  such  models  has  also  be  reported  in  the 
literature  18-20  Again  for  the  sake  of  completness  we  have  derived  the  model  kinetics  and 
statistics  based  on  the  units  in  Fig  2,  in  Appendix  B.  We  shall  use  the  results  derived  there 
to  compare  with  the  predictions  of  the  combined  model. 

Among  the  past  work  in  this  area,  Riccardi  and  Williams  (along  with  DusekpO, 
have  recognized  the  approximate  nature  of  their  earlier  model  18  and  the  concept  of 
fragments  has  been  used  to  highlight  this  difference. 

C.  POSTGEL  CALCULATIONS:  Having  concluded  that  the  random  model  is  an 
approximation  to  the  chemistry  (equations  3-5),  we  shall  extend  only  the  combined  or  the 
Bokare  and  Gandhi  model  to  predict  the  post-gel  properties  of  the  amine-epoxy  systems. 
For  the  postgel  calculations  the  set  of  superspecies  remains  the  same  as  in  the  pregel  stage 
and  the  same  kinetic  equations  are  used.  The  statistics  is  now  defined  in  terms  of 
expectation  values  of  discovering  a  finite  weight  attached  to  a  molecule,  when  a  unit  on  it, 
is  chosen  at  random. 

We  pick  an  epoxy  at  random  and  calculate  the  probability  that  the  weight  attached 
to  the  epoxy  looking  out  and  looking  in  is  finite  31. 


P(Fein)  =  [P(Feout)3(g-l) 


(24) 


Looking  in,  the  finite  probability  implies  that  the  weight  attached  looking  out  of  the  other 
'g-l)  epoxies  on  the  monomer  is  finite.  Looking  out  of  the  randomly  chosen  epoxy  the 
following  possibilities  exist: 

1)  It  is  attached  to  a  secondary  amine  sn  with  the  probability  :  p=nSn/e0  (Fig  la  group  E) 

2)  It  is  attached  to  a  tertiary  amine  tn  with  the  probability  :  p=ntn/eQ  (Fig  ib  group  E) 

3)  It  is  an  unreacted  epoxy  with  the  probability  :  p={  l-ZfaSn+ntnlJ/eo  =  e/eQ 

If  the  epoxy  chosen  at  random  is  reacted,  and  it  is  a  part  of  a  moelcule  with  finite  weight, 
the  weight  looking  into  all  the  units  attached  to  it  must  be  finite.  If  the  epoxy  is  unreacted 
the  probability  that  the  weight  looking  out  of  it,  is  finite,  is  unity. 


P(Ff ‘)  =  V  —  [{pCEin»n  l  *  Pdft]  +  V  *  PCPft] 

eQ  Lmi  e0 

n=l  n-1 


(25) 

Similarly  choosing  an  amine  at  random,  the  probability  that  the  weight  attachedto  it, 
looking  in,  is  finite,  if  the  weight  attached  looking  out  of  the  other  (f/2-1)  units  on  the 
monomer,  is  finite 


P(Fain)  =  {P(Faout)}(f/2-1)  (26) 

Looking  out  of  a  randomly  chosen  amine  we  again  distinguish  between : 

1)  secondary  amine  Sjj  ;probability  p  =  Sn/p0  (Fig  2a  group  A) 

2)  tertiary  amine  tn  ;  probability  p  =  tn/pQ  (Fig  2b  group  A) 

3)  unreacted  amine  p  ;  probability  p  =  p/p0  =  l-{Z(sn+tn)/p0} 


Then: 


p<C)  =£^{p<Ft">}D  *  + 

D»1  n*l 


i- 


/  -  M 

XV\ 

P„ 


Vn-l 


(27) 


In  the  first  term,  when  the  amine  chosen  happens  to  be  a  secondary  amine  sn,  the  weight 
looking  into  the  n  epoxies  on  that  chain  must  be  finite  and  similarly  if  the  amine  happens  to 
be  a  tertiary  amine.  Finally  the  last  term  reflects  the  fact  that  if  the  amine  were  unreacted  the 
probability  of  finite  weight  attached  to  it,  looking  out,  is  unity. 

These  are  a  set  of  four  non-linear  equations  24-27  in  four  variables,  P(Fein), 
P(Feout),  P(Fa*n),  P(Faout)  for  each  n.  We  again  use  probability  or  moment  generating 
functions  to  simplify  these  to  four  equations  for  all  n. 


s 


n  zn~'sr 


nz 


n-I 


(28  a-d) 


and  we  set  z  =  E(Fe'n)  in  28  a-d. 

We  have  already  obtained  s  and  t  as  a  function  of  z  from  our  kinetic  analysis  in  the 
pregel  stage 


(  X* 

s  _  zo  p  JL 
eo  X-l  LPo  *  lP°J  . 


(29) 


_t_ 

eo 


z2KAp 

(X-l)(l-y)(X-y) 


fp] 

r  p  / 

(1-Yl 

vPoJ  *(X“Y) 

lpo; 

J 

(30) 


Ka 

where  X=  -jt^I+KqZ) 


1  -z) 
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Performing  the  above  mentioned  steps  we  are  left  with  four  non-linear  equations 
24-27,  in  four  variables  that  can  be  solved  numerically  to  compute  the  finite  weight 
expectation  probabilities.  The  equations  and  the  numerical  method  are  delineated  in 
Appendix  A.  We  performed  the  numerical  analysis  to  predict  the  weight  fraction  of  sol  in 
the  post-gel  stage.  Other  properties  of  the  network,  e.g.  crosslink  density  can  be  similarly 
computed^  1-33  To  enable  the  reader  to  do  these  directly  we  also  present  the  finite 
expectation  values  for  two  interesting  cases. 


RESULTS 

A.  PREGEL  STAGE:  Since  we  wish  to  highlight  the  effect  of  etherification  reaction  in 
this  paper  we  will  ignore  the  substitution  effect  on  the  amine  groups.  There  are  two 
important  extremes  in  making  these  networks,  where  the  concerns  we  have  talked  about 
become  dominant 

1)  The  case  where  the  rate  of  etherification  is  faster  than  the  amine-epoxy  addition  reaction; 
that  is  k3>ki 

2)  The  case  where  the  network  is  formed  with  an  excess  of  epoxies  in  the  initial 
stoichiometric  recipe;  i.e.  r=2po/e0<l 

Under  both  these  operating  conditions  the  predictions  of  the  two  models  discussed 
would  be  significantly  different.  There  is  a  paucity  of  experimental  results  where  networks 
are  produced  under  either  of  the  above  mentioned  conditions.  There  is  a  commercial  epoxy 
resin  using  an  excess  of  epoxies  (DGEBA)  with  a  secondary  amine  (diethanonal  aminep4. 
We  are  currently  attempting  to  study  the  experimental  aspects  of  this  chemistry;  however, 
for  the  purposes  of  this  paper  we  shall  restrict  ourselves  to  the  theoretical  modelling 
considerations  of  such  networks. 

To  highlight  the  effect  of  etherification  rate  on  the  network  structure,  in  Fig  3  and  4 
the  weight  average  molecular  weight  profile  as  a  function  of  the  epoxy  conversion  is 
plotted,  for  both  the  models.  In  Fig  3  the  rate  of  etherification  is  ten  times  faster  than  the 
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amine-epoxy  addition  rate.  The  magnitude  of  the  deviations  between  the  two  models  is 
significant  In  Fig  4  where  the  etherification  rate  is  ten  times  slower,  there  is  still  an 
appreciable  difference  in  the  predictions  of  the  two  models,  but  as  the  relative  rate  of 
etherification  becomes  vanishingly  small  the  predictions  of  the  two  models  become 
identical.  In  both  Fig  3  and  4,  kj=k2  and  the  initial  stoichiometry  was  balanced,  that  is, 
2p0=eo  or  p=l/2.  Further  these  and  other  results  presented  here  are  for  bifunctional  epoxy 
(g=2)  and  bifunctional  amine  (f/2=2). 

In  Table  1,  the  gel  point  predictions  (as  the  conversion  of  the  hydrogen  in  the 
network  at  the  gel  point )  are  presented  for  the  two  models.  As  long  as  k3/k^<l  the 
differences  between  the  two  models  are  small,  but  as  the  rate  of  etherification  becomes 
significant  the  need  for  employing  the  correc :  combined  model  is  obvious.  In  Fig  5  the 
same  data  is  presented  in  a  graphical  form  for  interpolation  purposes. 

To  delineate  the  effect  of  imbalanced  stoichiometry  on  the  network  structure  in 
Table  2  and  3,  the  gel  point  is  presented  for  different  stoichiometric  ratios  r;  r  is  the  number 
of  hydrogens  to  the  number  of  epoxies  i.e.  2p0/e0.  In  Table  2  the  rate  of  etherification  is 
ten  times  slower  whereas  in  Table  3  it  is  ten  times  faster.  When  kj/Tq  is  less  than  unity 
and  r  is  greater  than  unity  the  predictions  of  the  two  models  are  within  a  few  percent  of 
each  other.  In  this  limit  the  epoxy-hydroxyl  reaction  is  almost  absent  and  the  predictions  of 
the  two  models  should  indeed  be  same.  Further  intuitively  one  would  expect  the  combined 
model  to  predict  a  faster  gel  point  and  a  sharper  growth  profile  as  it  faithfully  models  the 
chainwise  reaction.  This  too  is  observed. 

As  an  interesting  situation  consider  the  case  where  the  number  of  epoxies  are  four 
times  the  number  of  hydrogens  in  the  initial  recipe.  Further  assume  that  the  amine-epoxy 
system  is  such  that  the  rate  of  etherification  is  ten  times  slower  than  the  rate  of  amine-epoxy 
addition.  Then  Table  2  suggests  that  the  random  model  would  predict  that  even  at 
complete  hydrogen  conversion  the  system  does  not  gel,  whereas  the  combined  model 
predicts  a  gel  point  at  92%  hydrogen  conversion  and  27%  epoxy  conversion.  Under 
situations  like  this  where  the  rate  of  etherification  is  a  order  of  magnitute  less  than  the  main 
reaction,  a  random  model  would  still  be  inadequate  as  it  ignores  the  operative  growth 


chemistry.  Any  random  model  irrespective  of  the  statistics  it  employs,  either  combinatorial 
methods  probability  generation  functions  ^or  the  recursive  approach^  assumes  equal 
accessibility  of  all  ieactable  sites.  In  essence  it  makes  the  following  three  assumptions: 

1)  all  units  are  equally  accessible  to  the  reaction 

2)  all  units  combine  with  equal  randomness  irrespective  of  the  state  of  the  neighoring 
groups  or  the  size  of  the  growing  molecule 

3)  all  intramolecular  reactions  in  finite  species  are  ignored,  ie  no  loop  formation. 

In  chainwise  addition  reaction  where  the  reaction  is  conditional  upon  initiation, 
chains  do  not  grow  randomly  and  long  range  correlations  are  important  These  correlations 
change  the  nature  of  the  statistical  distribution  of  linear  chains  to  a  Poisson  distribution  13, 
assuming  all  initiations  occurred  at  the  same  time.  Assumptions  one  and  three  are 
unavoidable  at  present,  but  two  is  unrealistic  and  unnecessary.  It  is  for  this  reason  that  the 
need  to  understand  the  chemistry  of  the  network  growth  before  proceeding  to  build  an 
appropiate  model  for  it  is  important 

B.  POSTGEL  STAGE:  In  order  to  predict  the  various  network  structural  properties  in 
the  post  gel  stage,  the  probabilities  of  the  molecule  being  finite  .looking  out  of  a  randomly 
chosen  amine  or  epoxy  are  needed.  In  Fig  6  we  plot  these  finite  weight  probabilities 
looking  out  of  an  epoxy  for  the  cases  where  the  rate  of  etherification  is  ten  times  faster  and 
ten  times  slower.  In  Fig  7  we  plot  the  same  probabilities  for  the  amines.  In  the  case  of 
faster  etherification  rate  when  the  epoxy  conversion  is  99%  the  amine  conversion  is  only 
46%,  though  they  where  in  stoichiometric  balance  initially  and  the  finite  probability  looking 
out  an  amine  is  only  47%.  Such  systems  would  show  an  excess  of  unreacted  epoxy  even 
when  the  amines  are  driven  to  exhaustion.  Finally  in  Fig  8  the  effect  of  stoichiometric  ratio 
on  the  weight  fraction  solubles  during  the  postgel  stage  is  depicted.  The  weight  fraction 
solubles  are  given  by: 

ws  =  waP(Fa°ut)f/2  +  weP(Feout)g  (31) 
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where  wa  =  2pgMa/(  2pgMa+fMe)  is  the  weight  fraction  of  the  amines  in  the  network  and 
we  =  fMg/(  2pgMa+fMe)  is  the  weight  fraction  of  the  epoxies  in  the  network.  Again  the 
effect  of  imbalance  is  to  leave  an  excess  of  epoxies  unreacted  that  are  slowly  driven  to 
exhaustion. 

LIMITING  CASES  OF  THE  COMBINED  MODEL 

Since  it  contains  both  growth  mechanisms,  the  combined  model  can  be  used  to 
predict  the  network  structure,  in  either  of  the  two  operating  extremes.  If  the  network 
chemistry  is  purely  random  we  set  k3=0.  If  the  reaction  growth  is  entirely  anionic,  the 
amine-epoxy  reaction  can  be  considered  as  the  initiation  step  and  the  epoxy-epoxy  reaction 
as  the  propogation  step.  Further,  the  combined  model  overcomes  two  assumptions  often 
made  in  anionic  growth  models^: 

1)  the  initiation  step  is  instantaneous. 

2)  the  number  of  initiators  are  few  compared  to  the  monomers. 

In  the  combined  model  any  number  of  initiators  can  be  used.  Though  the  equations 
presented  here  are  limited  to  bifunctionally  paired  initiators  they  can  be  easily  extended  to 
the  general  case,  as  it  only  requires  modifying  the  recursive  analysis.  The  combined  model 
is  however,  limited  to  cases  similar  to  "living"  anionic  polymerization  since  it  contains  no 
termination  steps. 

SUMMARY 

We  have  compared  two  available  models  in  the  literature  for  treating  polymer 
networks  that  produce  a  secondary  site  which  reacts  through  a  different  growth  chemistry 
with  specific  application  to  amine-epoxy  systems.  We  have  explained  why  a  model  that 
operates  entirely  in  the  step  growth  paradigm  is  erroneous  and  that  the  exact  model  must 
consider  the  chainwise  nature  of  etherification  reaction  kineticaily.  We  have  cast  these 
differences  in  the  general  framework  of  superspecies.  We  have  shown  that  the  difference 
in  the  predictions  of  the  two  models  can  be  substantial  and  highlighted  such  situations. 
Further,  we  have  extended  the  combined  model  to  predict  the  network  properties  in  the 
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postgel  stage  of  the  amine-epoxy  systems. 
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APPENDIX  A  (POSTGEL  CALCULATIONS) 


First  we  give  equations  for  the  first  partial  derivatives  of  the  probability  generating 


functions  4a  &b. 


1  d(s+t)  -I 
e0  3lnz  J  z=l 
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+  w3 


(1-W2; 


:jpoj  -(wl-w2)|-  +  (wl-l)! 


where 


and 


'vf- 

A 


f+Ko 


2K, 


ka  1 

A  I 


(wl-1)  Kq(1-w2)  KQ(wl-w2)j 

(A2) 


wl  =  -A  (l-z+K^ 
2KA 

w2  =  — A(i-z) 

*0 


w3  = 


z2Kt 


(wl-l)(l-w2)(wl-w2) 


z  =  P(Wein) 


The  solution  to  equations  24-27  is  obtained  by  assuming  a  value  of  P(Fein)  which 
is  necessarily  between  zero  and  one.  Then  equation  27  can  be  solved  by  using  equations 
A1  and  A2  with  z=l  and  z=P(Fein)  for  the  last  and  the  first  two  terms  respectively.  26  can 
be  solved  using  27.  To  solve  25  the  first  partial  derivatives  of  the  probability  generating 
functions  are  needed  with  z=l  and  z=P(Fem)n*  *  which  are  given  by  equations  A1  and  A2. 
Having  solved  25,  24  can  be  solved  to  obtain  P(Fe*n).  This  value  is  compared  with  the 
assumed  root  and  the  analysis  repeated  till  desired  tolerance  is  achieved.  This  technique  is 
not  the  most  robust  mathematically  but  suffices  in  this  case. 
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APPENDIX  B  (  RANDOM  MODEL  ) 

Calculations  for  the  purely  Markovian  or  random  model: 

Notation: 

[N  i]  =  #  of  secondary  amines 

[N 12]  =  #  of  tertiary  amines 

[Epj]  =  #  of  epoxies  reacted  with  amine  only 

[Epl,2]=  #  of  epoxies  reacted  with  amine  and  epoxy 

[Ep2i2l=  #  of  epoxies  reacted  with  epoxy  and  epoxy 

[Ep2]  =  #  of  epoxies  reacted  with  one  epoxy 

[N0]  =  #  of  unreacted  amines  at  any  instant 

[Nfl  =  #  of  total  amines 

[Ep0]  =  #  of  unreacted  epoxies  at  any  instant 

[Epr]  =  #  of  total  epoxies 

f  =  #  of  hydrogens  on  the  amine  monomer 

g  =  functionality  of  the  epoxy  monomer 

p  =  [NT]/[Epx]  the  initial  recipe  ratio 

x  =  [N0]/[NT] 

a  =  k3/2ki 


REACTIONS 

ki 

No  +  EP0  - >  N1  +  EPl 

k2 

NI  +  EP0  - *  N 1 ,2  +  EPl 

k3 

EPl  +  EP0  - *  EPl,2  +  EP2 

K 

Ep2  +  EPo  - ^  EP2,2  +  EP2 

Assuming  kj=k2  and  defining  a  =  k3/2kj  the  differential  equations  for  the  above 
reactions  can  be  solved  in  terms  of  x  the  concentration  of  unreacted  amines  in  the  network. 
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The  differential  equations  for  this  growth  kinetics  are: 


•  =  2t,tN0]fEp0i 
dJN) 

-  -jp-  =  -2k,lNJtEp0]  +  kJN^IEpo) 
d[N.  ,1 

=  -kl"W 

d[Ep.) 

-  =  yEPiMEW  -  2kj[N0][Ep0]  -  kJNjlOEpJ 
d[Ep,  2] 

-  •  -k3[Ep,][Ep0] 
d[Hp-,] 

-  ~~=  -^[EpJfEpJ 
d[Ep, J 

--^p-=  -k3[Ep2][Ep0] 
dfEp  J 

*  =  2kj[N0][Ep0]  +  kJN^tEpJ  +  ^[Ep.ltEpJ  +  ^[Ep^tEpJ 

(B 1-8) 


Two  cases  in  terms  of  a  emerge: 
1)  a*.5 

then 

N, 

pi=n^='2x  +27x 
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EPl,2 


Pi, 2  “  R_  P2 


Ep 


2,2 


P2,2  ~  F_  _P 


-4/7)-  2a lnx  -2  (l+2a) 


Po  -  I2--  1  -  <P,+  P2  +  PU  +  P2,2)  =  1-  P  (2  <!'2a>  +  2“  tax  +  (2 -4a)  Vxj 


(B9-15) 


2)  when  a  =  .5 


pj  =  -2x  +  2/7 

pjj  =  x  -  2/x  +  1 

pj  =  -  p  (/7  lnx) 

p2  =  p  (/7  lnx  -  2/7  +  2) 

pj,2  =  p(/7  (lnx -2) +2) 

p2  2  =  p  ( -lnx  +  4/7  -  /7  lnx  -  4) 

p0  =  1  +  p  lnx 

(B 16-22) 

Armed  with  these  probabilities  one  uses  statistical  arguments  to  build  the  network 
structure.  Using  the  recursive  technique^  the  expectation  value  of  weight  looking  in  and 
out  of  a  randomly  chosen  epoxy  and  amine  is  determined.  In  order  to  do  that  the  following 
probabilities  arc  also  needed: 

Pis  =  eP°xy  reacted  with  secondary  amine  only  =  77 

plt  =  epoxy  reacted  with  tertiary  amine  only  =  Pt  -  pls 

Pi  2  s  =  epoxy  reacted  with  epoxy  and  secondary  amine  =  pj  2  77 

Pj  2t  =  epoxy  reacted  with  epoxy  and  tertiary  amine  =  Pi  2  '  Pi  2 

(B23-26) 
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Then 


=  (Pj+  P,  2) ECO  +  (plt  +  p12  +  p,  2t  +  p2  +  2p2  2)  EOV*enp) 
E(W14)  =  Mep  +  (g-l)E(Wj) 

E0O  =  (pI  +  2pn)E(WIeJ) 


E(wil)  =  Mam  +  (f/2-l)E(wS 


(B27-30) 


Solving  these  and  substituting  for  the  molecular  weight  one  gets  the  following 
criterion  for  the  gel  point: 

1  =T1  s(g*l)[(PI  +pl,2)(f/2‘  1)(Pl  +  pO)  +  (P!t  +  Pl^,  +  Pl,2  +  P2+2P2,2j 


and  M„,  = 


2pgMam 


w~  fMep  +  2pgMamE^ Wep)  +  Mep  +  2pgMam^v"an,/ 


EO^) 


(B31-32) 


where  E(Wep)=gE(W^)  +  M 


ep 


(B33) 


E(Wam)  =  (f/2)E(W^t)  +  Mara 


(B  34) 


E(Wo^  (P^Pt^CMam^Cf/MKp^lM.p)  4-  (Plt4'P1)2t+P2+Plt2+2P2.2)Mep  (B35) 


BOO  =  (pI+2pa)(Mep+(g-l)E(W^)) 


(B  36) 
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FIG  la:  General  structure  of  a  secondary  amine  ^ 

FIG  lb:  General  structure  of  a  tertiary  amine  tn 

FIG  2:  Smallest  statistical  units  needed  in  the  random  model  to  knit  the  network  fabric 

FIG  3:  Weight  average  molecular  weight  vs  epoxy  conversion  for  both  the  random  and  combined 
model;  balanced  stoichiometry,  r=l;  epoxy  functionality  g=2;  amine  functionality  f=4;  rate  of 
etherification  k3=  lk[and  k[=k2- 

FIG  4:  Weight  average  molecular  wight  vs  epoxy  conversion  for  both  the  random  and  combined 
model;  r=l,f=4,  g=2,  k3=10ki  andki=k2- 

FIG  5:  The  conversion  of  hydrogens  at  the  gel  point  vs  the  ratio  of  the  rate  of  etherification  k3  to  the 
primary  amine  -  epoxy  reaction  rate  kj;  r=l,  f=4,  g=2  and  kj=k2- 

FIG  6:  The  finite  probabilities  looking  out  of  an  epoxy  as  a  function  of  epoxy  conversion  for  two 
rate  ratios;  r=l,  f=4,  g=2  and  kj=k2. 

FIG  7:  The  finite  probabilities  looking  out  of  an  amine  as  a  function  of  epoxy  conversion  for  two 
rat^ ratios;  r=l,  f=4,  g=2  and  ki=k2- 

FIG  8:  The  weight  fraction  solubles  plotted  as  a  function  of  epoxy  conversion  for  different 
stoichiometric  ratios  r;  k3=10k3,  kj=k2,  f=4,  g=2. 
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Fig  2:  Smallest  statistical  units  needed  in  the  random  model  to  knit  the  network  fabric 
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TABLE  1 


Conversion  of  hydrogens  at  the  gel  point  for  different  relative  rates  of  etherification 

f=4,  g=2,  r=l. 


k3/k1 

Combined 

model 

Random  model 

%  difference 

.001 

.58 

.58 

~0% 

.01 

.56 

.56 

0s- 

O 

t 

.1 

.53 

.55 

3.8% 

1 

.39 

.47 

20% 

1  0 

.18 

.25 

39% 

100 

.05 

.09 

80% 

1000 

.01 

.03 

200% 

TABLE 2 


Conversion  of  amine  hydrogens  and  epoxies  at  the  gel  point  for  different  stoichiometric  ratios  r. 
The  rate  of  etherification  is  10  times  slower  than  the  amine-apoxy  addition  rate. 


r 

initial 

COMBINED 

MODEL 

RANDOM  MODEL 

stoichiometric 

ratio 

hydrogen 

conversion 

epoxy 

conversion 

hydrogen 

conversion 

epoxy 

conversion 

1/4 

.92 

.27 

no 

gelation 

no 

gelation 

1/2 

.72 

.39 

.78 

.42 

1 

.53 

.55 

.58 

.58 

2 

.38 

.78 

.81 

.81 

%  DIFFERENCE 


hydrogen  epoxy 

conversion  conversion 


TABLE  3 


Conversion  of  amine  hydrogens  and  epoxies  at  the  gel  point  for  different  stoichiometric  ratios  r. 
The  rate  of  etherification  is  10  times  faster  than  the  amine-apoxy  addition  rate.  f=4,  g=2. 


r 

initial 

COMBINED 

MODEL 

RANDOM  MODEL 

%  DIFFERENCE 

stoichiometric 

hydrogen 

epoxy 

hydrogen 

epoxy 

hydrogen 

epoxy 

ratio 

conversion 

conversion 

conversion 

conversion 

conversion 

conversion 

1/4 

.28 

.19 

.47 

.53 

68% 

180% 

1/2 

.23 

.26 

.34 

.56 

48% 

115% 

1 

.18 

.35 

.25 

.63 

39% 

80% 

2 

.13 

.47 

.18 

.74 

38% 

57% 

